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I .  INTRODUCTION 


A  number  of  interior  ballistic  models  have  been  developed  for  the  regenerative 
liquid  oropellant  gun,  e.g.  the  models  developed  by  Gough1’2  and  Coffee.3  These 
and  similar  models  are  capable  of  accurately  simulating  the  overall  performance 
of  the  regenerative  liquid  propellant  gun  (RLPG) .  However,  some  portions  of  the 
interior  ballistic  process  are  not  treated  in  detail  in  these  models.  In 
particular,  the  liquid  injection  process  is  usually  modeled  using  the  steady 
state  Bernoulli  equation  with  a  correlation  to  account  for  the  pressure  drop 
across  the  orifice. 

The  authors  have  previously  presented  a  model  of  the  liquid  injection  process 
in  a  regenerative  liquid  propellant  gun.A~6  This  model  was  developed  to  investigate 
the  unexpectedly  high  values  of  the  liquid  discharge  coefficient  observed  in 
regenerative  gun  firings.76  The  model  is  based  on  an  adaptation  of  the  Lagrange 
approximation  to  the  liquid  flow  in  the  reservoir  of  an  RLPG. 

In  this  paper,  a  modified  development  of  the  model  equations  is  presented. 
The  primary  modification  is  the  inclusion  of  the  full  pressure  distribution  over 
the  entire  reservoir  in  the  calculation  of  the  space  mean  pressure  ,  and  subsequently 
in  the  model  equations.  In  the  earlier  version  of  the  model, *  the  injection 
orifice  and  the  region  near  the  head  of  the  injection  piston  were  excluded  in 
the  calculation  of  the  space  mean  pressure.  The  velocity  distribution  in  the 
liquid  has  also  been  modified,  eliminating  another  approximation  used  in  the 
earlier  version  of  the  model. 

The  resulting  model  equations  are  then  applied  to  a  simplified  RLPG  geometry, 
in  order  to  evaluate  the  impact  of  utilizing  a  Lagrange  gradient  model  in  the 
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propellant  reservoir.  In  a  companion  report,9  the  model  is  applied  to  the 
simulation  of  the  RLPG  interior  ballistic  process,  and  comparisons  with  exper¬ 
imental  data  are  presented. 

II.  DEFINITION  OF  THE  CONTROL  VOLUME 

The  liquid  injection  model,  as  currently  formulated,  is  applicable  to  the 
regenerative  gun  configurations  known  as  Concept  VI  and  Concept  VIA,  shown  in 
Figures  1  and  2  respectively.  The  interior  ballistic  process  is  initiated  by 
firing  an  igniter  which  pressurizes  the  combustion  chamber.  The  chamber  pressure 
acting  on  the  injection  piston  forces  it  to  the  rear,  compressing  the  liquid  in 
the  reservoir.  After  an  initial  transient  period,  the  pressure  in  the  liquid 
reservoir  will  exceed  the  combustion  chamber  pressure  as  a  result  of  the  dif¬ 
ferential  area  across  the  injection  piston.  As  the  injection  piston  moves  to 
the  rear,  opening  the  injection  orifice,  liquid  propellant  will  be  injected  into 
the  combustion  chamber. 

Concepts  VI  and  VIA  are  very  similar  in  design.  The  primary  differences 
are  the  elimination  of  the  contour  at  the  rear  of  the  control  rod  in  Concept  VI 
(which  is  used  to  decelerate  the  piston  at  the  end  of  its  travel)  and  the 
introduction  of  a  damping  fluid  at  the  rear  of  the  injection  piston  shaft  to 
control  its  motion  in  Concept  VIA .  In  the  development  of  the  model,  we  will 
focus  on  Concept  VI,  but  the  simplifications  in  applying  the  model  to  Concept 
VIA  are  straightforward.3 

The  control  volume  for  Concept  VI  is  illustrated  in  Figure  3.  The  contours 
of  the  piston  and  the  reservoir  are  approximated  by  straight  line  segments  as 
indicated.  The  control  rod  and  transducer  block  are  fixed  in  the  reference  frame 
of  the  chamber.  The  origin  of  the  coordinate  system  fixed  in  the  chamber  frame 
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of  reference  is  ac  che  rear  (left  hand)  end  of  the  reservoir,  and  x  is  the 
coordinate  along  the  control  rod.  The  piston  moves  rearward  with  a  velocity 
u_(t),  and  the  points  s:(t),  s2( t >  ,  and  s3(t)  are  the  coordinates  of  the  fixed 
stations  on  the  inner  contour  of  the  piston  in  the  x  coordinated  system. 


BEFORE  FIRING 


CONTROL  INJECTION 

ROD  USTON 


Figure  1.  Concept  VI.  Regenerative  Liquid  Propellant  Gun 
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Consider  a  fixed  station  x  along  the  control  rod.  At  t  -  0,  the  corresponding 
point  on  the  inner  contour  of  the  injection  piston  is  defined  as  s(x,0).  At  the 
some  later  time  t,  the  piston  has  moved  to  the  rear  and  the  point  on  the  inner 
contour  of  the  injection  piston  at  the  fixed  station  x  is  now 

s(v.o  -  S(A'.O)  -  l‘up(i')dr  Cl) 

J  0 

such  that , 


<?s 

7t 


(2) 


and 


"  1  (3> 

Note  that  the  right  hand  face  of  the  control  volume  coincides  with  the  exit  plane 
of  the  injection  orifice,  s3(t),  such  that  the  control  volume  also  varies  with 
t  ime  . 

The  cross  sectional  area  of  the  control  volume  at  a  point  x  along  the  control 
rod  is 


d(v.f)  =  n[R\x.t)  -  r  ^  ( x )  ]  (4) 

where  R(x,t)  is  the  inner  radius  of  the  injection  piston  at  the  point  x  at  time 
t,  and  r „(.v)  is  the  radius  of  the  control  rod  at  the  point  x.  The  volume  of  the 
liquid  is  then  defined  by 
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-  ) 

1,(0=  |  A(x,t)dx  (5) 

J  C 

and  it  can  be  shown  that  the  rate  at  which  the  control  volume  changes  is 

0(0  -  -  «pM*  +  ^3(0]  -  -  U'AiiO  (6) 

where  AR  is  the  cross  sectional  area  of  the  reservoir  side  of  the  injection 
piston,  A3(t)  is  the  area  of  the  injection  orifice  at  the  point  s3(t)  ,  and 

■0(0=  -U+  d3(0. 


III.  LAGRANGE  APPROXIMATION  APPLIED  TO  THE  RESERVOIR 

The  equations  of  motion  for  the  liquid  in  the  reservoir,  continuity  and 
momentum  equations,  are  written  to  include  area  change  as  the  piston  moves 
rearward.  We  note  that  the  area  through  which  the  fluid  flows  is  a  function  of 
both  time  and  position,  since  the  contoured  piston  moves  rearward  over  a  contoured 
bolt.  The  one-dimensional  equations  of  motion  are  then 

<*(P  d)  +  <Hpi -A)  Q 

dt  +  dx  (  ' 


and , 


a(pt.i)  a(pt2d)  dp 

-  *  -  =  -  A  — 

d  I  dx  dx 

where  p.i.A  and  P  are  all  functions  of  both  position  and  time. 


(8) 


The  Lagrange  assumption,  density  is  a  function  of  time  only  and  is  thus 
constant  over  the  control  volume  such  that  the  spatial  derivative  is  zero,  is 
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a  good  approximation  in  the  case  of  the  LP  reservoir  since  the  liquid  density 
only  varies  by  about  4%  over  the  entire  ballistic  cycle  and  the  spatial  variation 
over  the  reservoir  at  any  given  time  is  much  less  than  this.  Therefore,  applying 
the  Lagrange  approximation,  Equation  (7)  becomes 

A  1  £?  =  .  iA  _  (9) 

pi?/  3  I  3  x 

Noting  that  p  =  mL/i  , ,  where  mL  is  the  mass  of  liquid  remaining  in  the  reservoir 
at  time  t  and 

mL  -  -  p  i3  ((3  +  up)  (10) 

is  the  mass  flux  out  of  the  control  volume,  and  using  Equation  (6)  we  obtain, 


1  3  p  _  ((3d  u  i>  d  k )  (11) 

p  dt  l  g 

The  numerator  on  the  right  hand  side  of  Equation  (11)  is  just  the  difference 
between  the  volume  of  liquid  exiting  the  reservoir  per  unit  time  and  the  volume 
swept  out  per  unit  time  by  the  piston  motion,  [(t3  +  up)A3  -  up(d3  +  d,)].  If 

the  liquid  were  assumed  to  be  incompressible,  these  terms  would  be  identical  by 
definition  and  dp/dl  =  0 

The  time  derivative  of  the  area,  A(x,t),  is  evaluated  using  Equations  (2) 
and  (4),  giving  us 


a  A  ( x- .  ( ) 
31 


dA 

3  A 

dS 

3t  i 

= 

(12) 


or , 


7 


it 


-  u , 


3[nR2(x.t)] 

dx 


(13) 


Using  Equations  (10)  and  (12)  in  Equation  (9)  we  obtain, 


dvA 

dx 


-  U  , 


ajn j?2(x.Q]  +  Jvla,  -  UPA, 
1  dx 


A(x.t) 


(14) 


and  upon  integrating  we  have 


t(-v,0  A(x.t)  =  up  {n[Rl  -  R\x,t )]}  ♦  (v2A3  -  u„  A  .  (15) 

v  i 

where  R0  is  the  inner  radius  of  the  injection  piston  at  the  origin,  i.e.  the 
rear  wall  of  the  reservoir. 


Rewriting  Equation  (15)  we  have 


t(.v,0  /I(x.f)  =  up  a(x.t)  *  (r3/i3  -  upAt ) 


1'U.Q 

V. 


(16) 


where , 


a(x.f)  »  -  n  Rz(x.t) 

such  that 


(17) 


A(x.t)  +  a(x.f)  -  [nl?o  -  nr*(x)]  .  (18) 

This  velocity  distribution,  Equation  (16),  exhibits  the  correct  behavior 
at  the  boundaries  of  the  control  volume, 

v(0 ,t)  -  0  (19) 
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and 


t(s3.f)  =  i. 


(20) 


where  we  have  used  the  fact  that, 


a(s3./) 


(21) 


Using  Equation  (7),  the  momentum  equation,  (8),  can  be  rewritten  in  the 


form 


dP 

d  X 


di 


di 


pi  7“  -  P  t: 

d  X  oi 


(22) 


Integrating  on  [0,x]  we  have 


P(x.l)  =  PD(t)  -  ^pi2(x,/)  -  p  £  u(x'.t)dx'  (23) 

which  is  a  form  of  the  unsteady  Bernoulli  equation.  Evaluating  Equation  (23) 
at  s3,  we  have 

P3(0  =  P  0(O  -  i  PJ3  -  P  f  v(x.t)dx  .  (24) 

C.  J  0 

also  a  form  of  the  unsteady  Bernoulli  equation,  where  P3(t)  is  defined  as  the 
pressure  at  the  right  hand  boundary  of  the  control  volume,  and  is  just  the  gas 
pressure  in  the  combustion  chamber. 


The  space  mean  pressure  in  the  liquid  reservoir  is  then  defined  by 


P{t)  A(x  ,l)dx  =  P(x,t)A(x,l)dx 

Jo 


p 

Jo 


(25) 


s 


or 


P(0 


_1_ 


x .  I )  A  ( x .  0  d  x 


Using  Equation  (23)  ,  we  obtain 


p0(o  -  no 


v2(x,  l)A(x.t)dx 

A(x .()  f  t i(x',t)dx'dx 
Jo 


where  P(t)  is  obtained  from  the  equation  of  state  for  liquid. 


(26) 


(27) 


Equation  (27)  is  utilized  in  Equation  (23)  to  determine  the  pressure 
distribution  in  the  reservoir,  and  in  Equation  (24)  to  complete  the  definition 
of  the  unsteady  Bernoulli  equation  for  the  reservoir.  We  now  consider  the 
integrals  involved  in  these  equations. 

1.  PRESSURE  DISTRIBUTION: 


In  order  to  evaluate  the  integral  in  Equation  (24)  we  first  differentiate 
Equation  (16)  to  obtain 


v(  \  .t)A(x.t) 


v(x.t)A(x.()  =  upa(x,t)  *  upa(x  ,t) 


(upAt 


v3A3 


Uj  A 3 ) 


V(x.O 

V, 


A  3) 


U(*.Q 

V, 


V*Y(x.Q 

V\ 


(28) 


where 
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A(x.t) 

da  (x.t) 

p  dx 

(29) 

a  (x.t) 

d a (x  A) 

=  Up  dx 

(30) 

A3(t)  = 

dnrl(x) 
u-  i* 

(31) 

l  '(•(,() 

-  ~  Up  a(x.t) 

(32) 

and 


t',(0  =  -  upAL(t)  .  (33) 

Substituting  Equations  (29)  through  (33)  in  Equation  (28)  and  rearranging  terms, 
we  obtain 


A,  j  Al  T  a(x,  t)  Al  K(x,<)~|\ 

Up[UpTL  t3^J\V^;[/l(x,£)_  ^(*-0  V,  J  j 


(34) 


Now,  using  Equation  (34),  the  integral  in  Equation  (23)  becomes 


j"  v(x'.t)dx'  -  Up 

-M  *  •  0  -  “  P 

r  ^ 
LUp^ 

3 

J2(x.l) 

*  Up  Z,(x.Q  + 

lip  J3(x,t)  - 

.  /*3' 

L2(x,t) 

+  upt3 

A 1 

+  L“p^i  " 

^3 

-/*(*.<) 

• 

(35) 

where , 


rx  a(x'.f)  da(x.t)  , 

— -  - - — ; — -dx 

0  42(x'.0  dx 


(36) 


J2(v./) 


J3(.v.O 


J4(.v.O 


Js(.V,0 


r  *'(** 

.0  7lz  ia(x',!)j_  , 

Jo  h,  A*(x‘,t)  dx' 

r  1 

da(x'.t)  ,  , 

—  - dx 

Jo  d  (  X  , 

0  ax' 

dnr£(x) 

/•*  l'(x'.l)  dx' 

dx 

Jo  r,  d(x',o 

al  r* 

a ( x ' .  0  >4  £  K  (  x ' ,  f) 

l’t  Jo 

_/!(*'. 0  4(x'.0  l/,  . 

(37) 


(38) 


(39) 


(40) 


£,(*.0 


r’  ^Lindx- 

0  A(x'.l ) 


(41) 


and 


M*-0 


^  f  ^iOdx- 

V',  Jo  d(x'.0 


(42) 
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Note  that  the  integrals  Jj  through  J5  are  dimensionless,  while  the  integrals  L: 
and  L2  have  units  of  length.  Rearranging  terms  in  Equation  (35)  and  substituting 
in  Equation  (23),  the  pressure  distribution  in  the  liquid  becomes. 


P(v.f)  -  MO 


-pt2(x-.0  -  pu2  [J,(x,0  +  J3(x,t)] 


V*AL 

PM'a  l-0(X-0]  -  pup 


[J2(x.O  -  J5(x,t)] 


0(*-0  -  -7-L2(.x't) 


pt 


A  3 
al 


(43) 


2.  UNSTEADY  BERNOULLI  EQUATION: 

Having  evaluated  the  integral  in  Equation  (23),  the  corresponding  integral 
in  Equation  (24)  can  be  determined  in  a  straight  forward  manner.  The  integrals 
defined  in  Equations  (36) -(42)  are  evaluated  at  s3,  and  are  written  in  the  form 

J,(s3.0  =  J°\t)  (44) 

and 

Z,(s3.0  =  L°3(t)  (45) 

etc.,  such  that  Equation  (24)  becomes 


p30)  -  P o(0  - 


2P*-3  -  P“p  (-O')  +  ^33(0] 


p“, 


p^pt's  t^53(0]  -  pu; 


vro)  -  JVto] 

L°\0  -  ^I“(0 


Pl'3 


A-.  m 

^f(o 


(46) 


3.  EVALUATION  OF  P0(t) : 


In  order  to  determine  P0(t) ,  which  is  required  to  complete  Equations  (43) 
and  (46) ,  we  must  evaluate  the  two  integrals  on  the  right  hand  side  of  Equation 
(27).  Using  Equation  (16),  the  first  integral  can  be  written  in  the  form, 


where 


--  /  «?(x.f).4(.v.l)d.v  =  ^  [■/?(<)] 


-  pup  UP—  -  t'3^  [J“3(0] 


P  f  ^  *  ^3  |  r  f  03  /  j  \  i 

2  p.i,  '  tJe  ( )] 


=  2 


(47) 


-On 


V, 


a\xA) 

d(x.n 


dx 


(48) 


o 


£i  f’3  i(v.Q  q(y.Q 

l7,  Jo  K,  d(x,t)  * 


(49) 


H 


Al 

A(x.t) 


dx 


(50) 


The  second  integral  in  Equation  (27)  is  easily  evaluated  by  noting  that  it  is 
simply  the  volume  average  of  Equation  (35)  over  the  reservoir.  Applying  this 
observation  to  the  integrals  in  Equations  (36) -(42)  and  defining 


and 


jV  (o 


i  rs- 

—  I  A(x  ,t)J  t(x  ,t)dx 

l  a  -  o 


(51) 


7?J(0  =  .—  /  A(x.t)L,(x.t)clx  (52) 

l  JP  Jo 

etc . ,  we  have 

I  Mx.t)  )*  i(x.t)dx'dx  -  pu2p[7f(T)  +  J§3(l)) 

l  /P  J  o  ^  0 

-  papL^  -  [7f(7j  -  7f(T)] 

L  al  al J 

_  r _  ^  _ ' 

+  pup«3[j;3(o]  *  p^P  *?3(0  -  -j-  i°3(n 

*  PI  3  X  W) 

_  ■'>  i 

-  p  i)7(0  •  (53) 

The  pressure  at  the  rear  wall  is  then  determined  by  substituting  Equations  (47) 
and  (53)  in  Equation  (27), 

Pjn  -  P{>)  +  !r*(0  *  pI7T(0  •  (54) 


15 


IV .  MOMENTUM  EQUATION 


The  momentum  equation  for  the  control  volume  shown  in  Figure  3,  in  the 
reference  frame  of  the  chamber  is, 


Ahu.  +  ~r  f  i’pdl  *  f  i’Pi  d  A  =  -  f  Pd  A  ,  (55) 

oi  J  c\  j  cs  j 

where  d A  is  the  outward  directed  normal  from  the  element  of  control  surface. 
Rewriting  Equation  (54)  we  have, 


a 

Jt 


\  -  2  - 

piAdxjt  *  pi-3d3i  = 


I  -s 

iJ. 

[PaA7  -  PcfA,  -  P3(A„  *  A,)].' 


(56) 


where  the  control  volume  has  been  extended  to  include  the  piston  shaft;  Pcf 
the  pressure  exerted  by  a  control  fluid  on  the  area  of  the  piston  shaft,  As;  Ap 
is  the  area  of  the  chamber  face  of  the  piston;  and  Ar  is  the  total  area  of  the 
transducer  block,  including  the  forward  axial  projected  area  of  the  damping  taper 
on  the  bolt  in  Concept  VI.  The  term  which  accounts  for  the  control  fluid  is 
used  for  the  Concept  VIA  configuration,  and  will  be  ignored  for  the  moment. 
Thus,  Equation  (55)  becomes, 


MPup 


7,  {/; 


P  V  Ad  X  )  —  P  3  (.  A  p  +  -^3)  ^0-^7  + 


(57) 


where  the  second  term  on  the  left-hand  side  is  just  ~{m/t;(l)} 
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We  will  first  evaluate  the  integral  in  the  above  equation  and  then  take  the 
time  derivative.  Using  Equation  (16),  we  have 


P 


r‘- 

j  i  (x . I) A( x.  I )d x 


a  (x .  t)dx 


-  P (Upl| 


l'3»3> 


1 


V* 


K(x ,  t)dx 


(58) 


or 


p  i(v.l)l(.v.l)d.v  -  mLur  J,(<) 


m  i(acdjf  1/3^3) 

— L°3\t )  (59) 


where 


=  —  )  ti  (  v  .  / )  d  v 

1  c 


(60) 


/  3  (0  -  r-  !  I  ( v./)d  v 


jp  -  0 


(61) 


Now,  taking  the  derivative,  we  have 


^  /  fs  \ 

^iPj  <(x.OA(x.Odxj  -  mtap 


d»3(0  -  77^(0 

1 


*  ml(3  -^£”(() 

L  ‘  * 


+PUP^3 


i”  (O 


<>nr*(v) 

ax 


♦  pup(tip^, +  ^3/t3) 


+  2piip(upd,  -  iV^-VtO 


j—  (“p  •l,-t3'43)[(Up'/1*-^3/,3)  +  UpAL)L°33(l)  . 

V  It 


(62) 


Substituting  in  Equation  (57),  we  obtain 
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U 


A 


*  t 


3,  rz°3(0 

1  t 


-  /»3(  -4  .^3)  -  P0Ar  +  pv23A3  -  pu  p(u  p  A  t  +  v3A  3) 


-  PU„«3 


r03/.^nr6J(r) 


LV  (O' 


+  2pup(iipd,-r343)./°3(Q 


*  upAL]L°3\t) 

1  * 


(63) 


which  is  the  force  balance  equation  for  the  piston 


V.  EQUATIONS  OF  MOTION 


The  unsteady  Bernoulli  equation.  Equation  (46),  and  the  force  balance 
equation  for  the  piston,  Equation  (63),  are  the  coupled,  ordinary  differantial 
equations  of  motion  governing  liquid  injection  and  injector  piston  motion.  These 
equations  can  be  rewritten  in  the  form, 


‘si 


«/. 


(0  - 


‘V 


*// 


(0 


1.-= 


-[P(t)-P3(t)]  -  -tj  -  */*(/) 


(64) 


and 


U'M'J'V)  -  t 3m;"(!)  -  ^(OM.+  Zi,]  -  /5(0dT-p^/)3-p[t/2d(0]  (65) 

respectively,  where  we  have  used  Equation  (54)  to  eliminate  P 0(<)  .  and  where 

-  i?3(0  (66) 

i',/(0  =  ~/.§3(0  (67) 
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f‘(Q 


u]\  *  fTFco- 


A, 


u,  l°l\t)--fL°23(t ) 
al 


T3i“CO 


(68) 


=  Al 


M  A(, 


j^i“3(0- j»3(0 

1  # 


(69) 


=  m  '  —1^(0 ! 


(70) 


f"  HO 


(0*2^6 *(0 


-  “p(“p^-t3^)>S3(0-J§3(fW?3(0] 


if  A* 


+  u„«3[js3(o]  -  1  [J° co] 


r- 


u„\LV( 0  -  -riS2CO 


.1 


,;L  ^ 


-^(ur^+l-3^3)  -  urv: 


i“J(0 


31  ^3(0 

anr^(x) 


ax 


p(u  pA  t  i3/A3)J9  (() 

u  p  ,4  ,  -  1 3  /!  3 


l 


[(up^-M3)  +  upal]L?0) 


(71) 


The  functions  t'2(f)  and  L'zA{t)  have  dimensions  of  velocity  squared  and  velocity 
squared  times  area  respectively,  as  suggested  by  the  function  names.  Note  that 
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these  functions  include  terms  involving  piston  and  liquid  acceleration,  up  and 
i3.  arising  from  the  evaluation  of  the  P0(O  in  terms  of  the  space  mean  pressure 

no. 


The  model  equations  have  been  incorporated  into  a  computer  model  and  applied 
to  Concept  VI  and  Concept  VIA  RLPG  configurations.  The  results  of  computer 
simulations,  and  comparisons  with  experimental  data  are  presented  in  the 
accompanying  paper.9 


VI.  DISCUSSION 


The  model  of  the  LP  reservoir  and  injection  process  which  has  been  developed 
in  the  preceding  sections  is  based  on  the  application  of  the  Lagrange  approximation 
to  a  flow  system  in  which  the  cross  sectional  area  varies  both  with  position 
along  the  direction  of  flow  and  with  time  at  a  fixed  position  in  the  system. 
As  a  result,  the  velocity  gradient,  Equation  (16),  is  quite  different  from  the 
more  usual  Lagrange  gradient  to  which  we  are  accustomed. 

Consider  the  development  of  the  standard  Lagrange  velocity  gradient  from 
the  continuity  equation, 


M-v-0 

dx 


(72) 


where  it  is  assumed  that  p  is  a  function  of  time  but  not  position  along  the  flow. 
Integrating,  and  assuming  that  v(0,t)-0.  we  have 


r(.v.<) 


(73) 


where 
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(74) 


"  _  P  1  =  l_  _  rn 

p  J  I'm' 

For  the  standard  Lagrange  gun,  in  which  the  propellant  mass  is  fixed  and  density 
varies  only  as  a  result  of  projectile  motion,  the  term  in  brackets  on  the 
right-hand  side  of  Equation  (73)  is  just  the  projectile  velocity  divided  by  the 
projectile  travel,  giving  us  the  familiar  form  of  the  linear  Lagrange  velocity 
distribution . 

If  we  now  permit  the  cross  sectional  area  to  vary  with  position  along  the 
flow  system,  but  assume  that  the  time  derivative  is  zero,  the  continuity  equation 

becomes 


dijx.t)A(x) 

6  X 

Integrating,  and  again  assuming  that  t( 0 . f )  —  0  ,  we  have 


P 

P  J 


A(x). 


(75) 


i(x.t)A(x) 


(76) 


or 


K-v.O 


_P  ’> '  ( -v ) 

p  J  A  ( x ) ' 


(77) 


From  Equation  (76),  we  can  see  that  the  mass  flux  at  a  position,  x,  in  the  flow 
is  proportional  to  the  time  rate  of  change  of  the  overall  system  density,  times 
that  portion  of  the  system  volume  from  the  origin  to  the  point  x  along  the  flow, 
i.e.  the  mass  flux  varies  linearly  with  volume  along  the  flow.  As  a  result,  the 
velocity  distribution  is  proportional  to  the  volume  to  area  ratio,  Equation  (77). 
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Finally,  we  permit  the  cross  sectional  area  to  vary  with  time  at  a  fixed 
position  in  the  flow,  as  well  as  with  position  along  the  flow.  The  continuity 
equation,  (9),  can  be  rewritten  to  give, 


<?t(.v,<)d(x.Q 

dx 


A(x.t)- 


<M(x,Q 


Integrating,  and  assuming  that  v(0 ,  t)  -  0  ,  we  have 


(78) 


t(x, Qd(x. 0 


V( 


ad(x'.t) 

dt 


dx' 


(79) 


The  first  term  on  the  light-hand  side  of  Equation  (79)  is  the  flux  arising  from 
changes  in  svstem  volume  (expansion  or  compression)  and  from  mass  loss  or  mass 
addition  to  the  system  as  a  whole.  As  noted  previously,  this  flux  varies  linearly 
with  volume  along  the  flow.  The  second  term  is  a  flux  arising  from  the  compression 
of  the  flow  due  to  a  change  in  flow  area.  If  the  integral  is  negative,  i.e.  if 
the  flow  area  between  the  origin  and  the  point  x  is  decreasing  with  time,  there 
is  a  positive  contribution  to  the  flux  at  the  point  x.  This  term  does  not  appear 
in  the  standard  Lagrange  gradient  model,  but  must  be  considered  here  due  to  the 
contour  of  the  moving  injection  piston. 


The  numerous  terms  which  appear  in  the  model  equations  are  all  traceable 
to  the  velocity  distribution,  Equation  (79)  .  The  role  of  the  velocity  distribution 
in  th*3  model  equations  is  summarized  in  Table  1.  The  terms  involving  v(x.t)  are 
related  to  the  kinetic  energy  of  the  fluid,  the  acceleration  of  the  flow,  and 
in  the  momentum  equation,  the  time  rate  of  change  of  the  average  fluid  momentum. 

Gough10  has  presented  a  development  of  the  ordinary  differential  equations 
of  motion  for  piston  motion  and  liquid  injection,  using  the  assumption  that  the 
liquid  is  incompressible.  The  resulting  equations  are  similar  in  form  to  Equations 
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Table  1.  Role  of  Pressure  Distribution  in  the  Model  Equations. 


Eq .  No . 

Equation 

Lq.  (2i) 

Pressure  Distribution 

jpe1(x.i).p/St/(A.\.')d.\' 

Eq.  (24) 

Unsteady  Bernoulli 

pfo3i>(x.l)dx 

Eq.  (27) 

Space  Mean  Pressure 

;pv2(t),piil(t) 

Eq.  (57) 

Momentum 

^{mlvCT)} 

(64)  and  (65)  ;  however,  Gough  has  approximated  the  integrals  appearing  in  Equations 
(24)  and  (57)  and  used  the  assumption  of  incompressibility  to  further  simplify 
the  solution. 

Comparing  Gough's  results  with  Equations  (64)  and  (65),  the  first  term  on 
the  left  hand  side  of  Equation  (64)  is  just  the  acceleration  of  the  liquid  in 
the  orifice.  Gough  has  approximated  this  term  by  v.  where  in  Gough's  notation, 
lH  is  the  length  of  the  orifice,  l  is  the  time  derivative  of  the  average  liquid 
velocity  in  the  orifice  and  is  "a  numerical  factor  close  to  unity".10  The 
function  i/ 2 ( f )  is  associated  with  the  kinetic  energy  of  the  liquid  approaching 
the  orifice.  The  term  in  Equation  (64)  involving  up  would  appear  to  be  related 
to  the  acceleration  of  the  liquid  approaching  the  orifice;  however,  there  is  no 
comparable  term  in  Gough's  development.10 

In  Equation  (65),  consists  of  the  piston  mass  and  the  mass  of  liquid 

which  is  being  accelerated  in  a  region  approaching  the  orifice.  Similarly, 
is  associated  with  the  mass  of  liquid  being  accelerated  through  the  orifice. 
Following  Gough's  arguments,  the  last  term  on  the  right-hand  side  of  Equation 

(65)  would  be  related  to  "a  correction  to  the  pressure  in  the  fuel  chamber  due 
to  the  velocity  of  the  fuel,"10  ie.  the  approach  velocity  of  the  LP. 
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VII.  SIMPLIFIED  RLPG  GEOMETRY 


The  coefficients  arising  in  the  model  equations,  Equations  (66) -(71),  are 
quite  complex,  and  it  is  instructive  to  consider  their  values  over  the  piston 
travel.  In  wider  to  tore  easily  evaluate  these  coefficients,  it  is  useful  to 
consider  a  simplified  RLPG  geometry.  Therefore,  we  will  consider  a  hypothetical 
variation  of  the  concepts  depicted  in  Figures  1  and  2,  in  which  the  bolt  has  a 
constant  radius  (i.e.  no  taper,  front  or  rear)  and  the  injection  piston  and 
transducer  block  are  configured  such  that  the  reservoir  and  injection  orifice 
are  right  circular  cylindrical  annuli.  The  functions  describing  the  area  and 
volume  of  the  reservoir  (analogous  to  Equations  (4),  (5)  and  (17)  )  are 


A  ( v 

.()  = 

diW(s,--v)  +  A3H (x  -  s2) 

(80) 

a(x 

■  0  - 

A,H(x-s2) 

(81) 

l'(v.O  -  ALxH(Sl-x)  *  [ALlt(t)+  A3(x-s2)]H(x-s2)  (82) 

where  /,(/)  is  the  distance  from  the  transducer  block  to  the  point  s,  on  the 
injection  piston. 

Using  Equations  (80) -(82),  the  integrals  required  to  evaluate  the  coefficients 
in  the  model  equations  are  easily  computed.  We  will  focus  on  the  effective 
lengths,  Equations  (66)  and  (67)  and  the  effective  masses,  Equations  (69)  and 
(70)  ,  since  they  depend  only  on  the  geometry  of  the  system  and  the  piston  position. 
The  integrals  required  for  this  evaluation  are, 
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(83) 


Table  2.  30-mm  Gun  Parameters 


1  j.  ( o ) 

7 

.4 

cm 

^2 

1 

.0 

cm 

al 

25. 

2 

cm2 

ar 

23. 

.3 

cm2 

A3 

1 

.9 

cm2 

Vr<0) 

188. 

.4 

cm3 

mL(0) 

269 

,8 

gm 

P  L 

1.432 

gm/ cm' 

Mp 

2109. 

.2 

gm 
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These  parameters  have  been  used  to  evaluate  Equations  (83) -(89),  which  have 
subsequently  been  used  to  calculate  the  effective  lengths  and  masses,  Equations 
(66),  (67),  (69)  and  (70),  also  as  a  function  of  piston  travel.  The  results  are 
presented  in  Figures  4  and  5 . 

The  behavior  of  £J//( <)  in  Figure  4  is  somewhat  surprising.  Its  initial  value 

is  about  3.325  cm,  or  about  45%  of  the  total  piston  travel.  It  decreases  steadily, 
becoming  negative  at  about  6.4  cm  of  piston  travel,  and  has  a  final  value  of 
about  -6.13  cm.  Thus,  the  effect  of  piston  acceleration  over  the  first  portion 
of  travel  is  to  increase  the  liquid  acceleration  t>3 .  while  the  effect  is  opposite 
over  the  last  1.0  cm  of  travel.  In  light  of  this  sign  change,  it  is  difficult 
to  attribute  physical  significance  to  this  quantity. 

The  second  effective  length  in  Figure  4,  L'J'(t),  has  an  initial  value  of 

about  1.27  cm,  or  about  1.27  times  the  orifice  length.  Its  final  value  is  0.50 
cm.  At  6.4  cm  travel,  the  value  of  is  still  approximately  1.0  cm.  The 
rapid  drop  to  the  final  value  occurs  over  the  last  1.0  cm  of  piston  travel.  The 
values  of  this  function  over  the  first  6.4  cm  of  travel  are  consistent  with  the 
approximation  used  by  Gough  in  reference  10. 

In  Figure  5,  we  see  that  the  initial  value  of  the  effective  piston  mass, 
M'/'( O-  is  about  5.8%  higher  than  the  actual  piston  mass.  This  represents  about 
122  g  of  liquid  which  is  being  accelerated  in  the  region  near  the  injection 
orifice,  or  about  45%  of  the  total  liquid  mass  at  the  beginning  of  the  injection 
process.  As  the  piston  moves  to  the  rear,  approaches  a  final  value 
slightly  less  than  1.0. 
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The  effective  liquid  mass,  in  Figure  5,  has  an  initial  value  of 

about  12.67  g,  or  about  4.66  times  the  mass  of  the  liquid  in  the  orifice.  As 
can  be  seen  in  Figure  5,  decreases  steadily  over  the  entire  piston 
travel,  reading  a  final  value  of  0.50. 

The  behaviors  of  these  four  coefficients  are  quite  reasonable  over  most  of 
the  piston  travel.  L'" (<)  and  do  begin  to  change  rapidly  very  close  to 
the  end  of  piston  stroke,  and  the  effective  length  coefficient,  L\n(t),  becomes 
negative  at  a  piston  travel  of  6.4  cm  when  the  ratio  of  the  reservoir  area,  AL  . 
to  the  reservoir  volume  begins  to  approach  1.  It  is  important  to  remember  that 
we  have  attempted  to  describe  a  flow  process  which  is  at  least  two-dimensional 
with  a  quasi  one  -  dimensional  model.  Over  most  of  the  piston  travel,  while  the 
reservoir  length  is  large  compared  to  its  radial  dimension,  this  approximation 
is  not  unreasonable.  However,  as  the  piston  approaches  the  rear  wall,  the 
two-dimensional  effects  become  dominant.  The  major  component  of  the  flow  velocity 
will  actually  be  in  the  radial  direction  near  the  end  of  piston  travel.  Therefore, 
we  anticipate  that  the  model  for  the  flow  in  the  LP  reservoir  and  injection 
orifice  which  has  been  developed  here,  might  not  be  valid  near  the  end  of  the 
piston  travel.  However,  as  will  be  demonstrated  in  the  following  paper,  the 
model  accurately  reproduces  the  experimental  LP  reservoir  pressure,  even  in  the 
latter  stages  of  piston  motion,  indicative  of  the  fact  that  the  correction  terms 
are  small  and  the  process  is  dominated  by  the  piston  inertia. 

VIII.  CONCLUSION 


A  model  of  the  LP  reservoir  and  LP  injection  in  a  regenerative  LP  gun  has 
been  presented.  This  model  is  a  revision  of  an  earlier  version,  the  primary 
difference  being  the  inclusion  of  the  full  pressure  distribution  in  the  model 
equations . 
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This  model  is  based  on  a  generalization  of  the  Lagrange  approximation  to 
address  the  variation  of  fluid  mass  in  the  reservoir  during  the  ballistic  cycle; 
the  variation  of  area  with  position  in  the  reservoir;  and  the  variation  of  area 
with  time  at  a  fixed  position  in  reservoir  due  to  the  rearward  motion  of  the 
contoured  injection  piston. 

In  a  companion  report,9  this  model  is  applied  to  simulate  the  injection 
process  of  two  RLPG  configurations  and  the  results  are  compared  with  experimental 
data . 
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A (x.t)  Cross  Sectional  Area  of  the  Flow 

a (  v  ,  ( )  nF  „  -  nF2( x ,  O 

■  hCO  d,+  d3(r) 

d;,  Cross  Sectional  Area  of  Chamber  Side  of  Piston 

A t  Cross  Sectional  Area  of  Reservoir  Side  of  Piston 

17  Cross  Sectional  Area  of  Transducer  Block 

d3(/)  Cross  Sectional  Area  of  Injection  Orifice 

//(v)  Heavyside  Function 

)  Non-Dimensional  Integral  Functions  Arising  in  Model  Equations 

J°’( f)  Space  Mean  Value  of  Jn(x,t ) 

/„(v .()  Function  with  Units  of  Length  Arising  in  Model  Equations 

1°3(()  Space  Mean  Value  of  I.(x'.f) 
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Liquid  Pressure  at  Transducer  Block 

Combustion  Chamber  Pressure 
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Radius  of  Inner  Surface  of  Piston 

R(O.t) 
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Point  on  Inner  Surface  of  Piston  at  Position  x  at  Time  t 
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